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1^^ ' Abstract 

O , . We examine the existence of relativistic stars in f{T) modified gravity 

^ ' and explicitly construct several classes of static perfect fiuid solutions. 

$_( I We derive the conservation equation from the complete f{T) gravity field 

Oij, equations and present the differences with its teleparallel counterpart. 

Firstly, we choose the tetrad field in the diagonal gauge and study the 

^vj , resulting field equations. Some exact solutions are explicitly constructed 

^ ■ and it is noted that these solutions have to give a constant torsion scalar. 

1/^ ' Next, we choose a non diagonal tetrad field which results in field equations 

1/^ , similar to those of general relativity. For specific models we are able 

^^ ' to construct exact solutions of these field equations. Among those new 

\l I classes of solutions, we find negative pressure solutions, and an interesting 

r~^ . class of polynomial solutions. 

o: 

T— I . 1 Introduction 

Modified theories of gravity have become very popular due to their ability to 

provide an alternative framework to understand dark energy. This is done by 

JH , modifying the gravitational Lagrangian to become an arbitrary function of its 

original argument, for instance f{R) instead of R in the Einstein- Hilbert action, 
see [D [2] • In General Relativity and its modifications one uses the metric 5^1/ 
and quantities derived from it to describe the gravitational field. 

There exists an equivalent formulation of General Relativity based on the 
idea of parallelism. Initially, Einstein aimed to unify electromagnetism and 
gravity based on the notion of absolute parallelism [3^, this however failed. 
Much later, the theory received more attention as an alternative theory of grav- 
ity which we now call the teleparallel equivalent of general relativity (TEGR), 
see [1]. 
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The idea behind this approach is to consider a more general manifold which 
contains in addition to curvature a quantity called torsion. The complete Rie- 
mann curvature tensor (part without torsion plus a contribution from torsion) is 
assumed to be zero and therefore one can in principle use either the torsion-free 
part or the torsion part to describe the gravitational field. The most convenient 
approach is to work with tetrad fields e*^ and a so-called Weitzenbock space, 
see [5]. The tetrad fields represent fields of ortho normal bases which belong 
to the tangent space of the manifold. This tangent space is Minkowski space 
equipped with the metric rjij and can be defined at any given point on the man- 
ifold. Note that e'^ has 16 components while the metric has only 10. However, 
the tetrads are invariant under local Lorentz rotations. 

Recently, modifications of TEGR have been studied mainly in the context 
of cosmology [6j. This theory is now known as /(T) gravity and is constructed 
with a generalised Lagrangian [7]. When compared with f{R) gravity, this 
modification is particularly appealing since its field equations are of second order 
and not of fourth order. Note that f{R) modified gravity can also be viewed 
as a second order system of equations when using the Palatini approach, see 
again [lj[2]- Most research on f{T) gravity is devoted to the theory's ability (or 
inability) to describe the observed accelerated expansion of the universe [S] , or 
comparing theory with observational data |10| . In addition, it has been shown 
that certain f(T) gravity models can provide us with a unification of early 
time inflation and late time accelerated expansion [TT] . Models which allow the 
equation of state to cross the phantom divide have also been found [12] . Other 
lines of research have also been followed, see for instance [13] . 

In two recent papers, static and spherically symmetric solutions were consid- 
ered [MllTS] in the context of /(T) gravity. In ^4J it was claimed that relativistic 
stars in f{T) do not exist, based on the general relativistic conservation equa- 
tion. In this paper we derive the conservation equation from first principles and 
show that it agrees with equation (18) of |14j. However, this equation should 
not be compared with its general relativistic analogue. We find solutions with 
constant torsion scalar T' = 0, similar to those found in [T3]. We proceed to 
consider various simple forms of ^^^ to solve the complete set of field equations, 
thereby showing the existence of relativistic stars. 

2 Teleparallel gravity and its modifications 

2.1 Basic equations and action 

The basic variables in the teleparallel approach to general relativity are the 
tetrad fields e*^ where the Greek indices (holonomic) denote the coordinates 
of the manifold while the Latin indices (anholonomic) denotes the frame. By 
staggering the frame and the coordinate index, we can use the same symbol for 
the matrix e*^ and its inverse. We define 

eVe.'' = S; , eVe/ = <5] . (2.1) 



We can define the metric via the tetrads by 

where rjij = diag(l, —1, —1, —1) is the tangent space metric, which is Minkowski 
space. Note that the determinant of the metric g is related to the determinant 
of the tetrad ^/—g = det(e*^) — e. The metric g is used to raise and lower 
coordinate indices and 77 raises and lowers frame indices. 

By assuming that the manifold is globally flat, the tetrad fields give rise to 
a connection defined by 

r;^ = e,^d,e\ = -eV9.e,^ , (2.3) 

which is the so-called Weitzenbock connection. Note that this connection is not 
the Levi-Civita connection since it is defined so that its torsion is zero. Since 
our manifold is flat, the notion of parallelism holds globally and therefore one 
speaks of absolute parallelism which is a synonym of teleparallelism. 
We define torsion and contortion by 

i^^^--^(r'^^-T'^^-T,'^"). (2.5) 

The contortion tensor can also be defined in terms of the Weitzenbock and Levi- 
Civita connections. It turns out to be useful to define the tensor Sa^"^ in the 
following way 

S-'^'^^l (^""^^ + SaTP% - SaT^'^p) ■ (2.6) 

Now we can define a torsion scalar T which is given by 

T = S^^-'T^, , (2.7) 

whose importance becomes clear in a moment. 

Due to the flatness of the manifold, we can express the Ricci scalar in the 
Einstein- Hilbert action in terms of the Weitzenbock connection, or equivalently 
torsion. This particular combination of torsion terms which appears in this 
context is the above mentioned T. It is thus rather natural to consider modi- 
fications of this action based on f{T) where / is an arbitrary function. Let us 
therefore consider the modified action (with geometrized units c — G — 1) 

5* = ^gravity + -Smattcr = -—- / 6 /(T) d'^X + / 6 Lmattor d'^ X . (2.8) 

2.2 Field equations and conservation equation 

Variations of the action (|2.8p with respect to the tetrads e'^ gives the field 
equations of f{T) modified gravity 

S.^^" hrd^^T -f e-^d^,{eSr)fT ~ T\,S,'''' Jt + Wf = in^" , (2.9) 



where Si'^'^ — ei'^Sa^" , Jt and Jtt denote the first and second derivatives of / 
with respect to T, see again [7j. Tfiu is the energy momentum tensor. In what 
follows we assume 7^^ to be an isotropic perfect fluid which is given by 

%u^pgtiu- {p + p)Uf_,u,j , (2-10) 

Conservation of the energy momentum tensor is ensured by the field equa- 
tions (I2.9P which we show explicitly in the following. Firstly, we rewrite the 
field equation in the form 

- eS;T\„S,''^^fT + ^S;/ = iTreS^T.-- . (2.11) 

We introduce the quantity ji" to which we refer to as a possible gauge current 
which represents the energy momentum of the gravitational field 

J^'' = -^ U''S,''-d^iT)fTT - efT\,S^''^fT + ie^v) • (2.12) 

If we set f{T) = T then ji^ reduces to the well known gauge current in telepar- 
allelism [TBI. Using ji^ , equation (I2.12p becomes 

e%a^(ee,"5,'''')/T - 4^ee%jV^ = ^^ed'^T.'' , (2.13) 

or equivalently 

d^iee^'^S^nfr - ^Trej," = ^Tree.PTp" . (2.14) 

Secondly, we now take the derivative of (|2.14p with respect to x"^ . The 
antisymmetry of S^^'^ (in the last pair of indices S'o-'"^ — —Sa'^^) implies the 
conservation equation 

47ra,(e(jT + m) = d^,{eSr)d,fT ■ (2.15) 

We immediately observe a few interesting facts. The presence of /(T) in the 
action affects the conservation equation in two different ways. On the one hand, 
the gauge current changes and on the other hand, also the right hand side is 
affected. When /(T) — T we recover the well-known conservation equation of 
TEGR [E] . In the following section we discuss static and spherically symmetric 
solutions in /(T) gravity and it will turn out that the contributions from the 
right-hand side term are crucial. 

3 Solutions with diagonal tetrad 

3.1 Field equations 

Consider the static spherically symmetric metric 

ds^ = e^'^^'^de - e''('')dr2 - R{rfdn^ , (3.1) 



where dfi^ = d6^ + sin^ 9 dip^ and where a, h and R are three unknown functions. 
One possible tetrad field (we can make arbitrary Lorentz transformations to the 
tetrads without changing the metric) can be written as 

eV = diag(e°('-)/2, e^('^)/2^ R{r), R{r) sin 9) , (3.2) 

to which we refer to as the diagonal gauge. The determinant of this diagonal 
matrix is the product of its elements so that we find e = e^'^+^^/^i?^ sin0. Based 
on this tetrad field, we can now write out explicitly the /(T) field equations (j2.9p 
in component form. To do so, we firstly compute the torsion scalar (J2.7I) and 
its derivative which become 

rv) ^ -2.- (-„":| + 4 _ ?^ , ^) _ r (. 4) , (.3) 

Inserting this and the components of the tensors S and T in to the diagonal 
components of the equation (|2.9I) yields 

(3.5) 



iirp = 




Anp — 


u^ ; 2 + 4 ' 


Airp = 


2 \2 Rj •'^^ 


- 


[T ^fR" a" fa' 
U+' U + 2+ 4 + 



In General Relativity the off diagonal components of the field equations vanish, 
but in /(T) gravity a diagonal tetrad field p.2p gives rise to an extra equatioi{3. 
This is the (r, 6) component: 

e^f COt(g)/TT ( ,i?'' , . ,,,i?'2 

„ii R I , I, I IIS. „R 

which can be written in terms of T' and gives 

e-^/2cot6' , , , 

-T'fTT = . (3.7) 



2i?2 

This implies that all solutions satisfy either Jtt = or T' = 0, where the former 
reduces the theory to TEGR. 



^ We are deeply indebted to Franco Fiorini for pointing this out to us and thereby correcting 
our previously incomplete field equations. 



The above field equations p.4p - p.7p provide us with four independent equa- 
tions for six unknown quantities, namely a{r), b{r), R{r), p{r), p{r) and /(T). 
Hence, this system of equations is under-determined, and in order to find solu- 
tions, we will need to make some reasonable assumptions to reduce the number 
of unknown functions to four. 

The most physical approach would be to prescribe an equation of state p = 
p(j>) which relates the energy density and the pressure, and to prescribe the 
function f{T). However, even in general relativity it turns out to be difficult 
to find explicit solutions for a realistic equation of state. Moreover, prescribing 
a possibly complicated function f(T) will make the field equations even harder 
to solve as they would contain more nonlinear terms. Thus, it turns out to be 
best to follow alternative routes. We will either make assumptions about the 
metric functions, or we will choose useful combinations of terms which simplify 
the field equations. 

3.2 Conservation equation II 

Based on the field equations p.4p - p.7p . we will re-derive the conservation di- 
rectly from them, without any reference to the previous results. In doing so, we 
will have two independent derivations of the conservation equation in the static 
and spherically symmetric setting based on a diagonal tetrad. Our approach 
here follows closely the well known derivation of the conservation equation in 
general relativity, see for instance [17] . In order to avoid working in TEGR, let 
us assume T' = in what follows. 



Differentiating p.Sp gives 



Wir) = -^fT. (3.8) 



An expression ioT p + p can be obtained using p.4p and (|3.5 



4n{p + p)=-^{R'{a' + b')-2R")fT, (3.9) 



while isotropy of the pressure implies 



f^-r\u^s:if^^^-^f^-^^ 



/y + :_ ( i:_ + r_+ ( r_ + i^^ (a'_6')^/j, = o. (3.10) 

\2R^ Aj-' 2 \R 2 \4: 2rJ^ ' ) '' ^ ' 

Let us multiply (|3.9p by a'/2 and (|3.10p by 2R' / R and subtract the two resulting 
equations. Since T' = we arrive at 

2na'{p + p) = ^jT. (3.11) 

Substituting (I3.1ip into (|3.8p leads to the /(T) conservation equation in a static 
and spherically symmetric spacetime with a diagonal tetrad field 

47rp' + 27ra'(/9 + p) = 0. (3.12) 



Notice that the right-hand side of this equation always vanishes because of the 
off-diagonal equation (|3.7I) which enforces JttT' = 0. Thus in /(T) gravity 
the static and spherically symmetric conservation equation in diagonal gauge 
coincides with its general relativistic counterpart, meaning that no physical 
differences can arise from it, compare with |14| . 

At this point we would like to recall the general conservation equation (|2.15p 
derived in the previous section 

-d^,{eSr)dAlT) + ^Trd^ieTn + ^^^,{ej^n = . 
Then, given the tetrad p.2p . wc find 

e'^/^ sin erf TT , a/2 2 ■ nf ,^^'P ^2P ^ ^'P 2p\ 

ine ' r amO [ p ^ \ 1 = . (3.13) 

2 \ 2 r 2 r J 

By rearranging and inserting T' = we see that this agrees with the conservation 
equation p.l2|) derived above. We are now ready to construct perfect fluid 
solutions for a constant torsion scalar explicitly. 

3.3 Solutions with T = 

Let us start by analysing the probably simplest solution. Let us assume that 
T = and insert this into (I3.3p . This fixes the function a{r) and we obtain 

a' = -^, a(r) = ln(ci/i?(r)) , (3.14) 

where ci is a constant of integration. Note that due to T = 0, / and its 
derivatives are now constants. Rewriting the field equations with this gives 

/ 1 e-^ fR' (R' ,, 



A'Kp = 



\2R^ 2 \R\R 

MO) /(o) 

2i?2 ^ 4 ' 
.-b /i?" R'2 y R' 



R"\\ 
2^ /t(0)- 


/(O) 
4 ' 


(3.15) 
(3.16) 






(3.17) 



while the conservation (|3.12p simplifies to 

4^P' + ^(P + P) = 0. (3.18) 

Isotropy of the pressure implies that h has to satisfy the differential equation 
e-^ / b' R' R" i?'2 \ 1 ^ ,^ „^ 



We can solve for h 



'-'"iTiJi^l- '^■^'" 



where C2 is another constant of integration. Therefore, the metric coefficients 
are given by 

and we arrive at the metric 

ds^ = 2ldt^ - ^^'^ dr^ - R^dn^ . (3.22) 

R C2-4:R ^ ' 

Having determined both metric functions, we can now attempt to solve the 
remaining equations. If we set friT = 0) = 0, then the field equations simply 
imply 

^0 = -Po = Y^ ' (3-23) 

thus we have found all unknown functions. Note that /(T) is arbitrary in the 
sense that only its value at the origin (and of its derivative) is of importance. 
Note that this excludes those /(T) which become singular as T -> 0. This 
solution yields a constant energy density po and pressure po, obeying the dark 
energy equation of state. Notice that the metric p.22p has a singularity when 
i? — > 0. Moreover, we expect a coordinate singularity at _R = C2/4 and it might 
be interesting to study this surface in more detail. 

3.4 Solutions with T' = 

Rather than assuming T = 0, we now assume T' ~ which is equivalent to 
assuming T = constant. These solutions are rather complicated but simplify 
substantially when considering the R{r) = r case. Let 

T{r) = To = constant , 

then the conservation equation p.l2p is again given by 

p' = -^ip + p). (3.24) 

The field equation p.Sp becomes 



To, ,^. fiTo) ^ MTo) 
y/TlTo) - -^- + ^^ 



47rp = —hiTo) — + -^-^ , (3.25) 



where fiTo) and /t(2o) are constants. As r — > 00 we have 



4^Poo = y/T(To) - ^ , (3.26) 



with Poo denoting the value of the pressure at infinity. 



Once we fix poo, p.26p can be viewed as a differential equations for /(Tq) 
wliose general solution is 

f{To) = fVn - 167rpoo , (3.27) 

where / is a constant of integration. Note that in this model Poo plays the role of 
a cosmological constant. However, though it represents a very simple spherically 
symmetric solution, it seems not to be of general interest in cosmology due to 
its incompatibility with standard teleparallel gravity. 
Using the equation for T, 

^0 = ^ («'^ + 1) ' (3.28) 

in order to write a' in terms of b and comparing p.4p and p.Sp we obtain the 
following differential equation for b, 

where Y{r) = e~^^'^\ An analytical solution to this differential equation seems 
hard (if not impossible) to find. However under certain approximations we can 
arrive to some solutions. 

First, consider the r^Tp <C 1 regime where p.29p reduces to 

Y' Y 2 

27 + 2^ + ^ = 0- (^-^0) 

The solution is then 

y(r) = e-^M = ^-4, (3.31) 

r 

which of course coincides with (|3.21l) . the solution for T — Q. To find also a{r) 
we can use (|3.28p . which in the regime r^T^ <C 1 gives again the solution p.2ip . 
Thus, at zeroth order in t^Tq the general solution matches the T = solution, 
as one would certainly expect. 

On the other hand, we can analyze the opposite regime: r^T^ 3> 1. Equa- 
tion p.29p becomes now 

^r' + ^y + r2To = 0, (3.32) 

and the solution is given by 

Y{r) = e^W = ^ " Y-^ (3-33) 

where fci is a constant. Again, to find a{r) we go back to p.28p which for 
r'^To » 1 reduces to 

«' - ^ - (3-34) 



«W = -oln ^ -]^ (3-35) 



and gives 



with fc2 another integration constant. Thus the metric in the regime r^To 3> 1 
is 

This metric is singular for r — )■ and r — >■ oo and presents an horizon for 

r=(2fci/ro)i/8. 

3.5 Triviality of the Einstein static universe 

Consider the case R{r) = r and the metric (13. ip with e"'^'"' and e'''^'"' fixed such 
that 

ds^ = dt^ V^tir^ - r^dn^ , (3.37) 

so that the corresponding diagonal tetrad field is given by 

e\ = diag (l, -7f=^^ r, r sin^") . (3.38) 

Since we have now chosen three functions, the system of equations is closed. 
For this choice the torsion scalar reads 

T = -a(l_t!). T--t, (3.39) 

and the field equations p.4p - p.6p become 



4(1 - kr^) ^__ , f^ ^\ ^_ , / 
2(1 - kr^) , f ^ A, f 



4^Po^ ^ :, V TT+(;^-fc)/T + T^ (3.40) 

^^P--(7^-k]fT~^, (3.41) 

4^P = - ^'^T' ' fTT-(T^-k]fT-^, (3.42) 



e-^/2cot6' , 
= ^^, T7tt . (3.43) 

The last field equation, and the isotropy of the pressure respectively imply 

T'fTT = , _ 2(l-fcr^) ^^^ ^ Q _ ^3 44^ 

Since 1 — kr^ cannot be zero, this can only be satisfied if frr = which takes 
us back to TEGR. Note that we cannot achieve T' = due to p.39p . 
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4 Solutions with off diagonal tetrad 

4.1 Field equations 

Consider again the static spherically symmetric metric 



ds^ 



.-{r)de - e'^'^Ur^ - Rirfdn^ 



(4.1) 



where d^^ — dO'^ + iiin 9 dip^ and where a, b and R are three unknown functions. 
Another possible tetrad field can be written as 



/ e°/2 \ 

e''/^ sin 6 cos 4> R cos 9 cos (p ~R sin 9 sin (f) 

e^/^ sin 6 sin i? cos 6 sin </) i? sin 6 cos 

V e''/2cos6' ~Rsm9 



(4.2) 



T^W 



("+'')/2i?(r)2sin6l. The 



see for instance |16j . The determinant of e*^ is e = e 
torsion scalar and its derivative are 



i?2 



(4.3) 



T'(r) 



,-6/2 



(4i?" + 2i?' (a' - 6') + i? (2a" - a'6')) 



+ ^ [RR"a' + i?'2 (a' - b') + R' {2R" + R {a" - a'b'))) - ^^ . (4.4) 
R'^ R 

Inserting this and the components of the tensors S and T into the equation (j2.9p 
yields 



A-Kp ■ 



-6/2 



R 



-{R'e 



^-b/2 



l)T'f'_ 



TT 



T 1 \ 
4^2^/ 



h 



,-b 



+ ^ {2RR" - RR'b' + R") /t - { , (4.5) 



4^P-(i-|-|^^'(^' + i?«'))/. + {, 



(4.6) 



47rp = — - 



R ' 



eV2 



^-i? + ki«" 



r /tt 

i?' Ra' , , , 
y + — )(«'-&') 



i?a' 



{• (4-7) 



The above field equations (|4.5p - (|4.7p give three independent equations for 
our six unknown quantities. As before, this system of equations is under- 
determined, we will reduce the number of unknown functions by assuming suit- 
able conditions. Note that there is no equation enforcing the constancy of the 
torsion scalar in this non diagonal gauge. 
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4.2 Conservation equation III 

We will derive the conservation equation one more time, now for the off diagonal 
tetrad. Taking the derivative of (|4.6p gives 

4^(0 = 1^ {2R' [e"/' Ra' R') + Ra'e''^) T' Jtt ^It 
- 1^ (f {R'{a' + h') 2R") + {a' {h'R' - R") - a"R') + ^) fr • (4.8) 
Next, we take a combination of the field equations (|4.5p and (|4.6p to obtain 

Mp + p) = ^ {R'e-'/^ - l) T'/tt - ^ (i?' + i?«') /t 

+ 2^ (2i?i?" - RR'h' + i?'2) /t . (4.9) 
Isotropy of pressure allows us to write 

We now multiply (|T^ by a72 and (HTTU|) by 2i?7-R and subtract the two re- 
sulting equations 

2^a\p + p) =-^ {2R' (e^/2 - Ra' R') + Ra'e"/') T' Jtt + ^/t 

comparing this with equation (|4.8p . we obtain the conservation equation for the 
off diagonal tetrad field (|4.2p 

p' + Y(p + p)=0. (4.12) 

It is surprising that one recovers the general relativistic conservation equation di- 
rectly in this approach. This seems to indicate that the non diagonal tetrad (|4.2p 
is a very good starting point to study static and spherically symmetric perfect 
fluid solutions in /(T) gravity. It also eliminate the problematic field equation 
in diagonal gauge which posed stringent constraints on the possible solutions. 

4.3 Solutions with fe = 

Consider our metric (|4.1I) and tetrad field (|4.2p with b{r) = and R{r) — r, 
then we have the metric 

ds^ = e'^'^'-Ut^ - dr^ - r^dn"^ , (4.13) 
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and for all values of a and R the torsion scalar is given by 

r(r)=0. (4.14) 

Again, since we have a vanishing torsion scalar, / and its derivatives are con- 
stant. The field equations simplify to 

4np=-l^, (4.15) 



4 
2r-' ' ' 4 



47rp = -^/T(0) + ^, (4.16) 



4-I' = -(^ + T(^ + -')yM0) + ^. (4.17) 



4 4 V 2 r 
Isotropy of the pressure yields 

«" + ^---0, (4.18) 

2 r 

which we can solve for a and find 

a = 21n(r2 + ci) + lnc2, (4.19) 

where C\ and ci are constants of integration. Thus we arrive at the metric 

ds^ = C2(r2 + c^fdt^ - dr^ - r^dfl^ . (4.20) 

Since / and fx are constants we can label /(O) = /i and /t(0) = f2- Using this 
and the metric coefficient a, we can write the field equations more explicitly 

(4.21) 

(4.22) 

(4.23) 





Anp 


h 

4 ' 








4:TTP 


2 
r2 + 


-/2 + 

Cl 


./l 

4 ' 




Anp 


2 

r2 + 


-,/2 + 
Cl 


h 
4 ■ 


lately, we see 


that p and 


p are given by 






Po = - 


h 
IGtt' 








P= - 


1 
27r(r2 + 


—J2 + 
Cl) 


A 

167r 



(4.24) 
(4.25) 
Notice that the pressure is regular everywhere provided that ci > 0. 

4.4 Solutions with /^(r = 0) = 

We now build a class of solutions for R{r) = r. Constraining the metric function 
b to be 

6 = 21n(l + ra') , (4.26) 
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leads to T = for all a. The field equations (|4.5l) - (|4.7p reduce to 

^ fT{0){2ra" + a'{ra'{ra' + 3) + A}) /(O) 
2r{ra' + lf 4 ' 

4,,.^^ + M. (4.2«) 

2 (r-^a' + r) 4 

_ /^(0)(a-(ra--l)(m- + 2)-2ra") /(O) 
8r(ra' + l)^ 4 ' 

If we further assume /t(0) = 0, so that we only consider particular /(T) models, 
these field equations imply 

The constraint /t(0) = is satisfied by a large amount of /(T) models. For 
example, we can take /(T) = T" (with n 7^ 1) or /(T) = cos{kT). The first 
gives p = p ^ 0, while the second gives p = —p = I/(167r). For all these models 
we have an infinite class of solutions given by choosing an arbitrary form for 
a{r) and constraining b{r) with the condition (J4.26I) . For instance, considering 
a Schwarzschild form for a{r), 

a(r) = In (1 - i^/r) , (4.31) 

we find 

b{r) = -2\n{l- K/r) , (4.32) 

with K being a constant. The metric becomes 

ds^ = (1 - K/r) df - (1 - K/ry^ dr^ - r'^d^^ . (4.33) 

which is very similar, but not equal, to the Schwarzschild solution. 

We can compare these solutions with the ones we found in Sec. (|3.3|) with 
a diagonal tetrad. We notice that assuming /t(0) = implies a Dark Energy 
equation of state in both cases. This suggests that, for these models, such a 
characteristic does not depend upon the choice of the tetrad. 

4.5 A negative pressure solution 

Following the route of Sec. (|4.4|) we now build a new interesting solution which 
applies to every /(T) model. We still assume R{r) = r and b{r) as in (|4.26p 
which again gives T = 0. As a particular form for a{r) we consider 

a(r) = --ln(fcr), (4.34) 

5 
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where fc is a constant. This implies b(r) — —2 In 5, corresponding to the foUowing 
metric, 

ds2 = (k r)-^l^dt^ - — dr2 - r'^dn'^ . (4.35) 

Note that the factor 1/25 in front of dr^ is physicahy unimportant. This is a 
special solution inasmuch as it renders both the field equations (j4.28p and (|4.29p 
equal to 

47rp = 1/(0) - ^/t(0) , (4.36) 

4 r^ 

without considering any constraint on f{T), which means that it works for all 
the /(T) models. The field equation (|4.27l) becomes 



1„,.. 12 
'4' 



47rp = --/(0) + -/T(0). (4.37) 



Note that if we require p to be positive everywhere we must have /(O) < and 
/t(0) > 0. This implies that p has to be negative for every r meaning that 
this is in general a negative pressure solution. Solutions of this type might be 
interesting in the context of dark energy or may have applications for wormhole 
like solutions. 

4.6 Solutions with T' = 

So far we found only solution having T = everywhere. In this section we show 
that a solution with constant T can be built for some specific /(T) models. 
Again we restrict the analysis to the R{r) — r case. Let us consider 

a{r) = \n{k r) ~ Tq r^ , 6(r)=ln4, (4.38) 

with k and Tq constants, corresponding to the metric 

ds^ ^kr e-^°''''df - Adr^ - r^dO^ . (4.39) 

Again we notice that the factor in front of dr^ is not of physical importance. 
The solution (|438l) immediately implies T = Tq. The field equations (|431) - (|4J| 
become 



>)+(t"8^^ 



inp = --/(To) + ( -^ - _ ) MTo) , (4.40) 



47rp = i/(ro) + ^/t(To) , (4.41) 



if (To) ^^ 



inp = -f{To) - ^ ,, ° ^ Mn) . (4.42) 



In order to satisfy the last two of these equations, we must require /t(2o) = 0. 
This is a constraint over all the possible /(T) models which implies that only 
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some specific models among these allow the solution (14.391) . From the field 
equations we obtain constant energy density and pressure, 

•^(-^") (A /IQ^ 

P=-P= -^— ■ (4-43) 

As an example of a f{T) model which allows for this solution we consider 
probably the simplest one: 

f{T)^T^^. (4.44) 

One can easily verify that /t(2o) = 0. For this model the energy density and 
pressure become p =^ —p = Tq/8. Note that requiring a positive energy density 
constrains Tq < 0, while in the approximation T <^ Tq standard teleparallel 
gravity (TEGR) is recovered. 

4.7 Polynomial Solutions 

Finally, we present a class of solutions which involves a non constant T and 
works for every model of the form f{T) = T^ /N. As before, we perform the 
analysis with R{r) = r. 
Consider the metric 

ds^ = {k r)"'dt^ - X^dr^ - r^dQ^ , (4.45) 

with n and m (real) numbers and k a constant of dimension 1 /length. The 
factor 1/n^ in front oi dr^ is physically unimportant. Metric (|4.45|) corresponds 
to choosing the two functions 

a{r) — mln{kr) , and &(?■)= 2 Inn, (4.46) 

which yield the non constant tensor scalar 

2(n - l)(-m + n ~ 1) 
T = -^ '-^-^^ '- . (4.47) 

Of course, we have to constrain n ^ 0, 7i ^ I and n ^ m + I in order to obtain 
a non-zero and regular T. Subtracting field equations (|4.28p and (|4.29p leads to 

{n~l){m~2n + 2){m-n + l) . {m^ - Am + An^ - 4} fr „ ,, ,„, 

n JTT H r^^ = ^ • y."^-^^) 

At this point we use expression (|4.47p in order to replace r^ with T in this last 
equation. In this manner we obtain a relatively simple differential equation for 
the function f{T) 

(m^ - 4to + 4n2 - 4) /t - 4(to - 2n + 2) T /tt = . (4.49) 



16 



If TO 7^ 2n + 2 the solution is 

/,T)..4r". wu„ ^^'^^^^. («o) 

and 7, ^ constants. Thus, for all these f{T) models metric (|4.45p represents 
a spherically symmetric solution. Note that TEGR is recovered when A^ = 1, 
which can happen for a wide choice of the parameters n and to,. For example, 
setting n — a/2 and to, = 2 gives N = 1. 

Finally, from the field equations we can read off the expressions for the 
energy density and pressure 

l&^p^~-f + £_hi{n,m)[^ '-^-^^ '-] , (4.51) 



Wnp = , + ^h,in,m) C^''-'^^7l^''-'^ Y , (4-52) 

where hi and /12 are two constants, depending on the parameters n and m in 
a complicated manner. If A^ > when r — >■ c» we have p = —p = j, while if 
TV < this happens at r = 0. The explicit forms of p and p do depend upon 
the choice of n and m. In general, we can state that a necessary condition to 
have both pressure and energy density positive everywhere is to set 7 = 0. 

Thus, we were able to build a class of solutions presenting a non-constant T, 
allowed by a large number of f{T) models of the form (|4.50p . These models are 
of general interest in cosmology and astrophysics since their study could lead to 
a better understanding of specific problems. 

5 Conclusions 

The objective of this paper was to show that relativistic stars exist in /(T) 
modified gravity and to derive some static and spherically symmetric perfect 
fluid solutions explicitly. Our starting point was a careful study of the field 
equations and the resulting conservation equations for two choices of the tetrad. 
We showed that the conservation equation of f{T) conceptually differs from 
those of General Relativity and its teleparallel equivalent. In doing so we also 
suggested a natural gauge current for f{T) modified gravity. We derive the 
energy-momentum conservation equation by using two different approaches and 
both these results agree. 

To begin with the study of exact solutions, we firstly studied some very 
simple models using a diagonal tetrad, by considering a constant torsion scalar 
T' = 0. Next, we examined some well known solutions of General Relativity, 
including Schwarzschild type solutions and the Einstein static universe. It turns 
out that the /(T) field equations for a Schwarzschild type ansatz of the metric 
are very difficult to solve and we were not able to make progress in this direction. 
We also found that the Einstein static universe with the usual assumptions 
of constant energy density and pressure only exists for the trivial f(T) = T, 
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which made it necessary to look for generahsations. A diagonal tetrad, though 
appealing, poses stringent constraints on possible solutions. 

We thus considered a non diagonal tetrad adapted to spherically symme- 
try For this choice of the tetrad, the conservation equation turns out to be 
equivalent to that of general relativity. We were able to construct a variety of 
exact solutions by making simplifying assumptions on the free functions. Among 
those solutions, we found one type of solution with a close resemblance to the 
Schwarzschild solution, we also found solutions with negative pressure. Finally, 
we derived new classes of solutions which we call polynomial solutions. This 
particular class of solutions has some very interesting properties. 
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